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Abstract This work aims to characterize mammographic x-ray beams incident and 
transmitted by breast phantoms (from 0 to 45 mm) composed from known proportion 
of glandular and adipose tissue-equivalent materials. This study was performed for 
mammographic x-ray beams generated by a mammography equipment using different 
target/filter combinations (Mo/Mo, Mo/Rh and W/Rh). It was studied the modification 
of spectra shape of the beams transmitted through different thicknesses of these 
materials. It was also evaluated the penetrability of these transmitted beams by its 
correlations to the HVL, which were experimentally estimated and derived from the 
x-ray spectra measured using a spectrometry system with a CdTe detector. The x-ray 
spectra transmitted by the phantom with higher density presented lower intensity than 
those transmitted by those with lower density, as expected. The differences between the 
HVL values derived from the spectra and those estimated using airkerma measurements 
are lesser than 6% for about 88% of the spectra measured in this work. The expected 
spectra variations with phantom thickness, revealed by the measured transmitted x-ray 
spectra, were also confirmed by HVL measurements and agree with the estimated 
attenuation curves. The motivation of the study was related to the robustness of the 
spectra as a descriptor of radiation beams and the possibility of using these transmitted 
spectra for dose assessment related to mammographic procedures. We can conclude 
that developed method is able to characterize mammographic x-ray beams making it 
possible the use of this kind of data for dose assessment in mammography. 
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Applications 

The x-ray spectra quantify the amount of photons and its energies in a x-ray Voi-5 No-i 

beam. It allows the calculation of derived quantities, generally applied for August2017 
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quality control and patient dose assessment, such as half value layer (HVL), 
air kerma, mean energy and others (Meyer et al. (2004)). An accurate 
representation of x-ray spectra used for radiological purposes is important 
for modeling the x-ray imaging process (Stumbo et al. (2004)). It allows the 
conversion of exposure information in dosimetric quantities for dose/image 
quality optimization (Bernhardt et al. (2006)). In the case of mammography, 
other applications of x-ray spectra have been reported in literature such as, 
radiation shielding calculation and evaluation of x-ray tube voltage (Kiinzel et 
al. (2004); Kiinzel et al. (2008)). Evaluation of mammographic radiation beam 
qualities and image properties have also been previously presented byTomal 
etal. (2013,2014). 

The knowledge of primary and transmitted x-ray spectra by materials 
simulators of breast tissues is very important once it can provide information 
about the composition and attenuation properties of these materials (Byng et 
al. (1998); Geraldelli et al. (2013); Johns and Yaffe (1987)). This information 
is essential for dose estimation. Moreover, these spectra can also provide 
information to characterize the source emitter of x-ray beams. 

Measurement of clinical mammography x-ray spectra present practical 
limitations mainly due to the high photon fluence produced by mammography 
equipments. Therefore, a common method to represent these spectra is using 
analytical or semi-empirical approaches (Blough et al. (1998); Boone et al. 
(1997); Hernandez and Boone (2014); Tucker et al. (1991)). However, this kind 
of spectrum becomes possible to be experimentally available after the advent 
of portable spectrometry system, which is capable to measure relatively high 
photon fluence. Mammographic x-ray spectra measured using these instruments 
generally requires complex alignment procedures, the use of pinholes and 
electrical or mechanical adjustments on the mammographic equipment. 
Different approaches have been reported in the literature (Matsumoto et al., 
2000; Miyajima and Imagawa, 2002; Stumbo et al., 2004; Abbene et al., 
2012; Zhang et al., 2012). Santos et al. (2017) have first introduced a study 
about the experimental x-ray spectra measurements in mammographic clinical 
equipments without changes in the original settings of these devices. Based on 
the Santos et al. (2017) methods the present work describes measurement of 
x-ray spectra to characterize mammographic x-ray beam transmitted by two 
different equivalent-tissue breast phantoms. The motivation of the study was 
related to the robustness of the spectra as a descriptor of radiation beams and 
the possibility of using these transmitted spectra for dose assessment related 
to mammographic procedures. Additionally, the possibility of measuring these 
spectra in clinical environment is relevant since it should represent realistically 
the radiation beams used in clinical practice. 
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This study was performed using mammographic x-ray beams generated 
using different target/filter combinations (Mo/Mo, Mo/Rh and W/Rh) incident 
on the breast phantoms composed by a known proportion of glandular and 
adipose tissue-equivalent materials. The x-ray spectra were measured using 
different thickness of the phantoms allowing the correlation of its morphology 
variation with the attenuator material. The results of this study can serve as a 
basis of future methods for dosimetry in mammography. 

2. MATERIAL AND METHODS 

Measurements of x -ray spectra and attenuation curves were performed in 
order to characterize the transmitted mammographic beams by different 
thicknesses of two tissue-equivalent breast phantoms. In this study, it was used 
x-ray beams generated using 28 kVp and each of the target/filter combinations 
(Mo/Mo, Mo/Rh and W/Rh) available in a Mammomat 3000 Nova (Siemens 
AG, Munich, DEU) mammographic equipment. Breast phantoms with two 
different tissue-equivalent compositions (CIRS Inc, model 012A) were 
positioned between the primary x-ray beam and the spectrometer, without the 
compression paddle, for measurements of transmitted spectra. These phantom 
compositions are equivalent to a breast 30/70 (30% of glandular and 70% of 
adipose tissue) and a 50/50 (50% of glandular and 50% of adipose tissue). In 
this work, the phantom equivalent of a 30/70 breast is referred as CIRS30/70 
and the phantom equivalent of a 50/50 is referred as CIRS50/50. For each 
tissue-equivalent composition, the slabs of the materials were superimposed 
intercepting the primary beam in order to assemble thicknesses between 5 
mm to 45 mm. Aluminum transmission curves were also measured for each 
target/filter and tissue-equivalent composition combinations. The flowchart 
presented in Figure 1 summarizes the tube parameters and experimental setup 
used for data acquisition. 
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Figure 1: Summary of the tube parameters and experimental setup used for the 
production of transmitted x-ray beam for which x-ray spectra and attenuation curves 
were acquired. 
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Figure 2: Spectrometer system and breast phantom positioned upon the breast support 
in a Mammomat 3000 Nova equipment, for transmitted x-ray spectra measurements. 

The transmitted x-ray spectra were measured using a spectrometer model 
XR-100T (Amptek, Inc., Bedford, MA) with a 9 mm 2 CdTe detector. The 
data processing was done by means of a digital pulse processor model PX4, 
also manufactured by Amptek. The detector was coupled to 25 pm to 100 pm 
collimator diameters, in order to reduce the detected photon flux and the scattered 
photon contribution on the detected beam. The collimator aperture was chosen 
depending on the thickness of the attenuation material. Figure 2 shows the 
spectrometer system positioned upon the breast support in a Mammomat 3000 
Nova equipment, for transmitted x-ray spectra measurement. At this position, 
the CdTe crystal was at 20.5 cm from the breast support. The lower base of 
the breast phantom was positioned above the spectrometer in a plastic support, 
such as a frame, at 32 cm from the breast support. For this mammography 
equipment, the source to the positioning assembly is fixed at 65 cm. 

The energy calibration of the spectrometer was done using standard 
radiation sources of 24l Am, 152 Eu and 133 Ba. The measured transmitted x-ray 
spectra were corrected by the detector efficiency, estimated using Monte Carlo 
method (Tomal et al. (2015)), Compton distortion and fluorescent escape 
fraction were reduced or removed using stripping procedure (Di Castro et 
al. (1984); Santos et al. (2016)). In this way, the corrected x-ray spectra are 
presented in units of fluence normalized by the current-time product (photons/ 
mm 2 .mAs.keV). 

For each x-ray spectrum, corresponding air kerma measurements were 
performed using a dedicated mammography ionization chamber model 10x5- 
6M (Radcal Corp., Monrovia, USA). This same detector was used to measure 
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Figure 3: Experimental setup for air kerma measurements for the derivation of 
attenuation curves and HVL values highlighting the chamber positioning in the light 
field (a) and in the x-ray field (b) of the mammography equipment, (c), presents the 
positioning of al foils and the breast phantom. 


transmission curves by aluminum. For these measurements, the x-ray field was 
collimated using a lead plate with a hole of approximately 1.6 cm positioned 
next to the tube collimators. Figure 3 presents the ionization chamber 
positioned for the air kerma measurements used to derive the attenuation 
curves and HVL. Figure 3 (a) presents the ionization chamber positioned at 4.5 
cm above the breast support in the light field of the mammography equipment 
delimited by the lead plate used. Figure 3 (b) presents the radiographic image 
of this detector positioned over the breast support such as shown in Figure 3 
(a). Figure 3 (c) presents the positioning of the aluminum filters and the breast 
phantom for air kerma measurements. 

The experimental data for the beam attenuation by aluminum were fitted 
using Archer’s model (Archer et al. (1983)) in order to estimate the transmission 
curves. The HVL values corresponding to each x-ray spectra measured were 
derived from these transmission curves. HVL values were also derived from 
experimental spectra and compared to those experimentally obtained using 
ionization chamber in order to assess the feasibility of obtaining such quantity 
from x-ray spectra. For this calculation, first the air kerma spectra were 
computed from the photon fluence spectra, N(E), as showed in Equation (1). 


N k {E) = N{E)x 


M„(£) 


xE 


(1) 
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Where is the photon fluence spectrum as a function of the energy, E, and are 
the mass energy transfer coefficients for air. Second, the air kerma transmitted 
by different thicknesses,*, of aluminum was computed analytically. The sum 
of these distributions is equivalent to the value of air kerma,, of a x-ray beam 
transmitted for different thicknesses of aluminum: 

K ar {x) = Y.N k (E)e 1 p -A E (2) 


Where are the mass attenuation coefficients for aluminum. Finally, the ratio 
was computed for each * value used. Then, the HVL values were obtained 
using cubic spline interpolation of values in order to find the * value (HVL) 
which satisfies the following relation: 


K ar {x = HVL ) 

*-(o) 


(3) 


The variation of x-ray spectra and its respective HVL values with the 
phantom thickness was analyzed for each anode/filter combination in order to 
characterize the typical beams transmitted for material that represents different 
thicknesses of breast tissue. 



Figure 4: HVL values obtained for transmitted x-ray spectra measured for each Mo/Mo, 
Mo/Rh and W/Rh target/filter combination as a function of the phantoms thicknesses. 
The filled and unfilled symbols represents the HVL values obtained for the x-ray 
spectra transmitted for the breast phantoms Cirs30/70 and Cirs50/50, respectively. 
These values represent the experimental HVL, estimated from measurements using 
ionization chamber. The maximum uncertainty obtained for these values was 5%. 
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3. RESULTS AND DISCUSSIONS 


The HVL values obtained from transmitted x-ray spectra measured for Mo/ 
Mo, Mo/Rh and W/Rh target/filter combination are illustrated in Figure 4. 
These results are presented as function of the phantoms thicknesses. The filled 
and unfilled symbols represents the HVL values obtained for the x-ray spectra 
transmitted for the breast phantoms CIRS30/70 and CIRS50/50, respectively. It 
can be noted that, for the same target/filter combination, the use of the phantom 
with higher glandular tissue-equivalent content intercepting the radiation 
beam results in transmitted beams with higher HVL values and consequently 
higher effective energy relative to those transmitted by the material with 
lower glandular tissue-equivalent content. It can also be observed that the 
transmitted W/Rh beams are the most penetrating and Mo/Mo beams are the 
less penetrating. The knowledge of the power of penetration of these beams is 
very important for consideration of energy dependence in dosimeters used for 
depth dose calculations. 

Figure 5 presents the transmission curves that represent the normalized air 
kerma for x-ray spectra transmitted through 5 to 45 mm of tissue equivalent 
phantoms CIRS30/70 and CIRS50/50 relative to the air kerma produced by 
the incident photons on the surface of these materials. These results were 
obtained for incident beams generated using 28 kVp different anode/filters 
combinations. This figure shows that the air kerma reduction occurs with 
increasing phantom thickness and this reduction is more pronounced for 
phantom with larger proportion of glandular tissue (CIRS50/50). The W/Rh 
radiation beams proved to be harder than Mo/Mo and Mo/Rh beams. These 
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Figure 5: Relative air kerma transmission by 0 to 45 mm of tissue equivalent phantoms 
Cirs30/70 and Cirs50/50. It was used incident beams generated using 28 kVp and Mo/ 
Mo, Mo/Rh and W/Rh anode/filter combinations. 
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Figure 6: Transmitted photon fiuence (photons/mm 2 .mAs.keV) by 0 to 45 mm of 
tissue-equivalent phantom Cirs50/50. The x-ray incident beam was generated using 
28 kVp and W/Rh target/filter combination. 

results are in accordance with those found with measured transmitted x-ray 
spectra presented in Figure 4. 

Figure 6 presents the photon fiuence spectra transmitted by 0 to 45 mm 
of tissue-equivalent phantom CIRS50/50 for W/Rh combination at 28 kVp. 
The spectra intensities presented in this figure decrease with increasing of 
material thickness. In this way, the spectrum with higher intensity corresponds 
to the primary beam, and that with the lower intensity corresponds to the beam 
transmitted by 45 mm of the attenuator material. Similar results were obtained 
for Mo/Mo and Mo/Rh target/filter combination for both breast phantoms used 
(CIRS50/50 and CIRS30/70). 

Figure 7 presents the comparisons between fiuence spectra transmitted 
by the same thickness of the breast phantoms used in this study. These 
comparisons refer to the transmitted spectra by 20 mm of the breast phantoms 
(Figure 7 (a)), 30 mm (Figure 7 (b)), 40 mm (Figure 7 (c)) and 45 mm (Figure 
7 (d)). The spectra transmitted by the material with composition equivalent 
to a breast 30/70 were plotted using black lines and those transmitted by the 
material equivalent to a breast 50/50 were represented using red lines lines. 

Table 1 presents the HVL values (mmAl) calculated from air kerma 
measurements performed using ionization chamber (Experimental F1VL) and 
those derived from the measured x-ray spectra. This quantity was calculated 
for all spectra measured in this study considering different target/filter 
combination and the several thicknesses of the breast phantoms used. The 
F1VL values derived from the spectra presented a difference to those estimated 
using air kerma measurements less than 6% for about 88% of the spectra used 
in this work. 
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Table 1: HVL values (mmAl) calculated from air kerma measurements performed 
using ionization chamber (Experimental HVL) and those derived from the measured 
x-ray spectra. 

CIRS30/70 

„ ^ , ... HVL Derived from spectra 

Experimental HVL (mmAl) (mmAl) 


Phantom 


thickness 

Mo/Mo 

Mo/Rh 

W/Rh 

Mo/Mo 

Mo/Rh 

W/Rh 

(mm) 







0 

0.312 (4) 

0.358 (2) 

0.481(9) 

0.326(2) 

0.371(2) 

0.498(2) 

5 

0.378(1) 

0.421(1) 

0.535(1) 

0.391(3) 

0.433(3) 

0.546(3) 

10 

0.418(1) 

0.460(9) 

0.567(3) 

0.421(2) 

0.462(3) 

0.578(3) 

20 

0.462(2) 

0.518(20) 

0.622(1) 

0.470(4) 

0.522(4) 

0.664(3) 

30 

0.500(4) 

0.552(2) 

0.673(4) 

0.501(4) 

0.568(4) 

0.674(4) 

40 

0.533(2) 

0.591(9) 

0.714(7) 

0.538(3) 

0.607(4) 

0.716(5) 

45 

0.545(1) 

0.612(6) 

0.734(7) 

0.549(5) 

0.621 (4) 

0.754(4) 

CIRS50/50 

0 

0.312(4) 

0.358(2) 

0.481(23) 

0.326(2) 

0.371(2) 

0.498(2) 

5 

0.373(3) 

0.417(2) 

0.513(4) 

0.393(3) 

0.436(3) 

0.542(3) 

10 

0.419(7) 

0.461(2) 

0.571(1) 

0.425(3) 

0.470 (3) 

0.570(4) 

20 

0.472(7) 

0.529(1) 

0.643(1) 

0.474(5) 

0.517(5) 

0.629(5) 

30 

0.524(17) 

0.577(5) 

0.695(30) 

0.495(8) 

0.539(7) 

0.658 (9) 

40 

0.555(25) 

0.619(4) 

0.736(2) 

0.501(11) 

0.536(10) 

0.669(11) 

45 

0.571(3) 

0.640(2) 

0.752(8) 

0.567(5) 

0.643(4) 

0.731(5) 
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Once the glandular tissue is denser than the adipose tissue (Hammerstein 
et al. (1979)) the spectra transmitted by the phantom with composition 50/50 
has lower intensity than those transmitted by the phantom with composition 
70/30, as expected. It can be noticed that the difference between the 
intensities of these spectra for the materials attenuator used increases with 
increasing phantom thickness. It occurs because of the ratio of the number 
of photons transmitted by the phantom CIRS 30/70 to those transmitted by 
the phantom CIRS50/50 for a given energy, should increase exponentially as 
a first approximation, in complying with the exponential attenuation law for 
monoenergetic beams. 
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Figure 7: Comparisons between x-ray spectra transmitted by some thicknesses of two 
different breast phantoms: (a) 20 mm, (b) 30 mm, (c) 40 mm and (d) 45 mm. 


CONCLUSIONS 

We tested a portable spectrometry system with CdTe detector to measure the 
x-ray spectra from radiation beams generated by a mammography equipment 
and transmitted by breast phantoms in order to characterize these spectra that 
represents the distribution of the unattenuated photons when breast of different 
thicknesses and composition were irradiated for examination. 

HVL measurement is a common method for characterization of x-ray beams 
for quality control and for assessment of equipment performance. In this study 
we have introduced a more sophisticated method to do this characterization, 
since from the x-ray spectra, many quantities of dosimetric interest, such as 
dose, HVL, mean energy, etc, can be derived. It is important to notice that 
the modifications expected to the transmitted beam for different target/filter 
combination and different attenuator material extracted from the measured 
transmitted x-ray spectra were consistently confirmed by HVL measurements 
and agreed with the estimated transmission curves. Therefore, it can be 


100 
























concluded that the developed method was able to adequately characterize 
mammographic x-ray beams. 

Once a method for transmitted spectra measurements is established, depth 
dose variation in these materials can be evaluated. Furthermore, these x-ray 
spectra can be used to derive absorbed dose and mean glandular dose at various 
depth of these phantoms. The knowledge of HVL and the effective energy 
transmitted by different tissue compositions and combinations of various 
target/filter also allow the development and characterization of appropriate 
detectors which facilitates more efficient production of clinical images. 
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